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Abstract. It is a common assumption amongst astronomers that, in the de- 
termination of the distances of remote sources from their apparent brightness, 
the cumulative gravitational lensing due to the matter in all the galaxies is the 
same, on average, as if the matter were uniformly distributed throughout the 
cosmos. The validity of this assumption is considered here by way of general 
Newtonian perturbations of Friedman-Robertson- Walker (FRW) cosmologies. 
The analysis is carried out in synchronous gauge, with particular attention to 
an additional gauge condition that must be imposed. The mean correction 
to the apparent magnitude-redshift relation is obtained for an arbitrary mean 
density perturbation. In the case of a zero mean density perturbation, when 
the intergalactic matter has a dust equation of state, then there is indeed a zero 
mean first order correction to the apparent magnitude-redshift relation for all 
redshifts. Point particle and Swiss cheese models are considered as particular 
cases. 



1. Introduction 

The purpose of this paper is threefold: first to give a general analysis of syn- 
chronous gauge Newtonian perturbations of Friedman-Robertson- Walker (FRW) 
cosmologies; second to obtain a condition which identifies a suitable homogenised 
FRW cosmology for comparison purposes, and third to obtain the relevant formula 
for the correction to the apparent magnitude-redshift relation for weak gravitational 
lensing. 

The classic paper of Lifschitz & Khalatnikov (1963) generated much interest 
in the study of perturbations of FRW cosmologies. Among the influential papers 
that followed were those of Hawking (1966) and Sachs & Wolfe (1967). A review 
that covers both classical and quantum aspects is Mukhanov ct al. (1992). One 
of the main topics in the literature is the growth of small perturbations in the 
early universe and the formation of the galaxies. The present work, by contrast, 
is concerned with the description of the universe as it is today, with the galaxies 
providing an essentially time independent perturbation of an FRW background. 

Although the use of a gauge invariant approach is frequently advocated, since 
it eliminates any need to specify and interpret gauge conditions, it is nonetheless 
more appropriate here to employ a synchronous gauge. This gauge is known to 
involve a certain ambiguity, as was pointed out by Mukhanov et al. (1992), and has 
given rise to difficulties of interpretation (Press & Vishniac 1980). It will therefore 
be necessary in the present work to identify an appropriate supplementary gauge 
condition which admits a natural interpretation. This will be done in §o. 
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Gravitational lensing has been considered by many authors and with different 
techniques (see e.g. Schneider et al. 1992; Bertotti 1966; Dyer & Roeder 1972,1973, 
1974; Kantowski 1998). One such technique has been the use of stochastic numerical 
integration of the lensing equations (Holz & Wald 1998; Dyer & Oattes 1988). The 
effects of weak gravitational lensing on the cosmic microwave background radiation 
have been considered by Seljak (1996) and by Sachs & Wolfe (1967). The optical 
properties of the Swiss cheese model of Einstein & Strauss (1945,1946) have been 
considered by Dyer & Roeder (1974). Both the Swiss cheese model and the point 
particle model of Newman & McVittie (1982) will be considered here as applications 
of a general theory. 

Newtonian perturbations of FRW cosmologies have previously been considered 
by McVittie (1931), Harrison (1967) and Newman & McVittie (1982). A central 
difficulty is that if one neglects terms involving the spatial curvature constant k, 
the total gravitational potential for a uniform distribution of particles will be every- 
where infinite. What does not seem to be widely appreciated is that this problem 
disappears when proper account is taken of either a positive or negative k. This 
will be discussed in §|| b- 

Much of the analysis will be carried out in the linear approximation in the grav- 
itational coupling constant K = 8irG/c 2 . Consequently only weak gravitational 
lensing will be described. In particular, shear is neglected because it has a second 
order effect on luminosity. Caustics are also neglected. These may both be sig- 
nificant restrictions. There are now well-known examples of astronomical images 
that show the effects of caustics, although there are as yet no firm estimates of the 
proportion of images for which caustics are involved. It has been claimed by Ellis 
et al. (1998) that light rays followed back from our present location in time and 
space meet a galaxy within a redshift < z < 5. If this is correct then strong 
gravitational lensing is likely to be significant for images with z > 5. 

The question as to whether, in general, weak gravitational lensing gives the 
same apparent magnitude-redshift relation on average as a best-fit FRW model was 
considered by Weinberg (1976) who gave two supporting arguments. The first was 
based on an analysis of double image lensing but was only valid for q <C 1, i.e. for 
<C 1. The second, subsequently generalised by Peacock (1986), was essentially 
based on photon conservation (Schneider et al. 1992, pp. 99-100), and was valid 
for all f2. Although widely accepted in the literature, these arguments have been 
criticised by Ellis et al. (1988), the first on the grounds that generic lensing gives 
rise to three images, not two, and the second at least in part because it effectively 
assumes the result to be proved. The criticism of Weinberg's first argument is 
presumably based on the odd number theorem of Burke (1981) and McKenzie 
(1985), according to which a transparent gravitational lens can give rise only to an 
odd number of images. However this theorem is flawed since, as was pointed out 
by Gottlieb (1994), it attempts to use 3-dimensional topology in a 4-dimensional 
setting. Nonetheless the criticism of Weinberg's second argument appears to stand. 
The present paper therefore gives a new supporting argument that, in the weak field 
limit, and when averaged over large angular scales, the apparent magnitude-redshift 
relation does indeed agree with that of a best-fit FRW model. 
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2. General perturbations of FRW cosmologies 
The FRW metric, in standard form, is 

g= -dt 2 + S 2 (t)h (2.1) 

O 

where h is a t- independent 3-metric given by 

Here and henceforth, indices i, j ■ ■ ■ run from 1 to 3, whilst indices a,b . . . will run 
from to 3. For the most part it will be possible to describe the level surfaces 
of t with respect to an arbitrary coordinate system x = {x 1 , x 2 , x 3 }. The space- 
time may then be conveniently described with respect to the coordinate system 
{x a : a = 0, . . . , 3} with x° = t. 

O O 

Quantities associated with g will carry a superscript . Covariant differentiation 

° o o o 

with respect to g will be denoted by ;, e.g. g a ^° c = 0. The energy tensor of g is 

T ab = Cp+ p)u a u b + °pg ab (2.3) 

where g ab is the inverse of g ab , u a := -~g ab t°i, is the unit future-directed normal to 
the level surfaces of t, 

*> = s|) + HH (2 ' 4) 

give the density p(t) and pressure p(t) respectively, and A is the cosmological con- 
stant, which is not assumed to vanish. A superscript dot denotes differentiation 
with respect to t. 

The Hubble and deceleration parameters for g are defined by 

^):=|| (2.6) 



S(t)S(t ) 
S 2 (t) 

respectively whilst the dimcnsionless density parameter is defined by 



(2-7) 



m 4^ • (2.8) 

3H 2 {t) 



Note that by ( |2.4[ ) and ( |2.5| ) one has 

jH{t) = -H 2 (t)(l+q(t)) = g£fi- 1<P + P) ■ (2-9) 

o 

Covariant differentiation with respect to h in the level surfaces of t will be de- 
noted by a subscript ||, e.g. hij\\k = 0. Quantities associated with the 3-metric h, 
apart from h itself will, in addition to the superscript , carry a prefix . Thus 
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the Laplacian associated with h will be denoted by ' 3 'A. The Riemann and Ricci 

o 

tensors and Ricci scalar of h are 

= 2k6 i [k h l]j (2.10) 
{3) Rij = 2khij (2.11) 
(3) R = 6k (2.12) 

respectively. 

The calculations of the remainder of this section pertain to a perturbation of the 



FRW metric ( |2.1|) of the particular form 

g = -dt 2 + S 2 (t)h (2.13) 

where 

h = h + 5h (2.14) 

is a 3-metric intrinsic to the level surfaces of t, i.e. transverse to u, and where Sh 
may depend upon all four coordinates. 

The scale factor S(t), the cosmological constant A and the spatial curvature 
constant k will all be taken to be the same as for the unperturbed metric g. Note 

o 

that since h is independent of t one has dthij — dtShij. Covariant differentiation 
with respect to h in the level surfaces of t will be denoted by ||, e.g. hij\\k = 0. 
Quantities associated with h, other then h itself, will carry a prefix ^ 3 - > . Thus 
the Laplacian associated with h will be denoted by and the Ricci tensor and 
Ricci scalar of h will be denoted by ^Rij and ^ 3 'R respectively. Indices i,j... on 
quantities associated with h will be raised and lowered by h _1 and h respectively; 
indices a,b . . . on quantities associated with g will be raised and lowered by g _1 

o 

and g respectively. Differences between quantities associated with g and g, or h 

and h, will be denoted by S, e.g. (5 (3) i?y := (3) R i3 - (3) i?y . 

The vector field u a = —t' a is irrotational and geodesic. The integral curves of u 
are to be regarded as the world lines of a preferred family of observers, described 
henceforth as comoving. In particular the galaxies are presumed to be comoving. 
The coordinate t thus carries physical significance. 

The second fundamental form of the level surfaces of t is given by 

Xab ■= u a . b . (2-15) 

Clearly \ab is transverse to u a and so defines a 3-tensor intrinsic to the level surfaces 
oft: 

Xa = 4" = ldt(S 2 (t)h ZJ ) = S{t)S(jt)hij + ^S 2 {t)d t 5h %3 . (2.16) 
The Hubble and deceleration parameters of u a with respect to g a b are defined by 
H := u a , a = i X a a = H{t) + U-> ;),<)!,,, (2.17) 

g:=ft(^)-l. (2.18) 
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Note that H and q are functions of all four coordinates. By (2.16) the Gauss 
equation for ^ R l jki leads to 



(3) i? = 6fc - S 2 {t)n{°p -p) + 2S 2 (t)( K ST tt - ^-h ij d t Sh i:j ) 



S(t) 



+ -ti^h^ n {d t 5h tj )(d t 5h mn ) . 



From (fyJI), {fiJ% , (pTlib and ( ^19| ) one has 



(2.19) 



«X« - ^d t (S 2 (t)Sh l3 ) 



5H = -V'dtShij 
SH 



H 



H 



H 



(2.20) 
(2.21) 
(2.22) 



6&R = 2S 2 (t)( K ST tt - ^-h^dtShij) + \h^ m h^ n {d t 5h t] )(d t 5h mn ) . (2.23) 

The components of the Ricci tensor of g ab are given by 
S(t) 



Rt 



Ru — Rti — h? ' {dt8hk{i)\\j\ 



S[l) 2 SHt) dt{S2{t)hl3dt5hl]) ~ l^iPt&hiM&tShu) (2.24) 

(2.25) 



1 A 



R tJ = (2k + 2S z {t) + S(t)S(t) + -S(t)S{t)h km d t 5h km )h l3 
+ 6 ^Rij - 2kSh ij + -L-atiSPtfiBtShij) 

+ ^s 2 (t)h km (d t 5h km )(d t sh t] ) - ^s 2 (t)h km (dSh lk )(d t 8h jm ) , (2.26) 

the last of these having been obtained by means of ( [2.11 ) . One also has 



R = 6l 



S 2 (t) S 2 (t) S(t)> S 2 (t) 



8&R 



(2.27) 



by means of (2.11) and 

S &R = h ij (3) i?y - ti ] (z) °R l3 = ti J {5 (3) % - 2k5h ij ) 
Hence the energy tensor T of g, as defined by Einstein's equations 



(2.28) 



nT ab = R ab - -Rg a b + ^9ab 



(2.29) 



CUMULATIVE GRAVITATIONAL LENSING 

has components 



KTtt = K ° P+ (3)R + Wj hij9t6hij + \( hijd t 6h ij) 2 ( 2 - 3 °) 

kT u = kT u = (hi k d t 5h k[i ) m (2.31) 



KTij = {S 2 (t) K p - -L-^d t (S 3 (t)h km d t 5h km ) 

1 



S 2 {t)hP^ m (d t 8h km ){d t 8h m )}h i3 



4 

+ *< 8 >iZ« - 2kSh ij - h^S^R+^dtiS^mo,,,, 



5hj 



+ 1 S 2 (t)h km {(d t Sh km )(d t 5h lJ ) - 2(d t Sh ik )(d t 5h jm )} (2.32) 



1 

t 

The Weyl tensor C a bcd of g a b may be decomposed in the standard form 

Cabcd = 8u[ a Eb][dU c ] + 2g a [ c E d ] b - 2g b [ c E d ] a - 2u e r/ eabf U[ c B d ] f - 2u e i] ecdf U[ a B b } f 

(2.33) 

where rjabcd is the alternating tensor, 

E ab := C acbd u c u d (2.34) 
Bab := *C acM w c U d (2.35) 

are the electric and magnetic parts of C a bcd, and 

*C a bcd ■= -jllabefC ef cd — -T] c defC ef a b (2.36) 

is the dual of C a bc d- 

In order to compute E ab one may use the identity 

2u b u a .[ bc ] = R dabc u b u d (2.37) 

to obtain 

E H = ~ \s 2 {t)(R tt + l -R)h l3 - ^ Xij - ^XikX k j (2.38) 
which by ( f2~16| ), dJI), ( pT26|) and ( |2~27| ) gives 
l S (t){d t (S(t)d t 6h i3 )-± 

^(S^R ij -2kSh ij -^h ij S ( - 3) R) 

3 

To compute B ab one may proceed directly from ( 2.3 5| ) to obtain 

cd 1 f 

Bab = -U e Ve C a(Xbd;c ~ -g c bRjdU } ) (2.40) 

which by the symmetry of B a b gives 

B ab = -u e ri e cd {a Xb)d,c ■ (2.41) 
By means of (2.16) one thus has 

Ba = -\s(t){ i3) v lm i(d t Sh jm ) lu + (3) v lm AdtSh im )^} (2.42) 



Eij = --s(t){d t (s(t)d t Sh tJ ) - -h l3 h km d t (s(t)d t Sh km )} 



+ -S 2 (t){h km (d t Sh km )(S t Sh lj ) - -h n vh km {d t 5h km ){d t 5h np )h v] } . (2.39) 
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where 

(3 W = ^u a Valjk (2.43) 

is the alternating tensor for hij . 

The calculations so far have all been exact. However, in order to facilitate further 
progress, the perturbation Sh will be treated as a power series in k that vanishes 
to zeroth order, and the perturbed metric g will be studied only to first order in 
«. Note that the spatial curvature constant k will not be treated as a power series 
in ft so, for example, terms with the coefficient kn will not be disregarded in first 
order approximations. This contrasts, in particular, with the work of Holz & Wald 
(1998). 

3. REDSHIFT AND EMISSION TIME 

Suppose a comoving observer with world line {x = Xi} in the perturbed space- 
time makes an observation at time t — t\ of a comoving source with redshift z\. The 
image seen by the observer is formed by a congruence of null geodesies centred on a 
null geodesic 7 connecting a point (to, xo) on the source world line {x = xo} to the 
observation point (ii,xi) on the observer's world line {x = xi}. For comparison, 
consider an observer in the unperturbed space-time who, at the same observation 
point (ii,xi), makes an observation of a comoving source with the same redshift 

00 
zi. In this case let 7 be the null geodesic, {x = x } the source world line and t 

o 

the time of emission. For definiteness let 7 be chosen such that its spatial direction 
at (ti,xi) coincides with that of 7. 

o 

The immediate problem is to compute the perturbation dto — to — to in the time 
o 
of emission. For g the emission time to is determined implicitly by 

l + (3-D 
S(t ) 

which for nearby sources gives 

ti-t = + 0(4) fOT 21 < 1 ■ (3-2) 

H(t x ) 

For g one may determine to as follows. Since 7 is a null geodesic one has 

dx l dx^ 
~dt~d~t 

d 2 x a „„ dx b dx c 



52 (*)*« — — = 1 ( 3 - 3 ) 



dX* =° ^ 

where A is an affine parameter along 7 with value Aq at (to,Xo) and value Ai at 



(ii,xi). By means of (3.3) the t component of (3.4) gives 



where 



d, f dt\ S(t) , . 

s 1 "U) + sM +KM ' ) = " (3 ' 5) 



,(„ := IsHmSH,/^ (3.6) 
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is defined along 7. The redshift as observed at time t is given by the standard 
formula 



l + z(t) = 



which combines with (gj) to give 

d 



dt 



In(l+z(t)) 



(dt\ \ 
(MA ) 

S(t) 



S(t) 



KV 1 if) 



Integration of ( |3.8j ) from to to t% yields 



1 + Zi = Wo) BM 



r (t)dt) 



t=t 



(3.7) 



(3.8) 



(3.9) 



which may be regarded as an implicit equation for to in terms of z\ and t\. To first 
order in k equation (3.9) gives 



St 



1 



H(t ) 



nv°(t)dt + 0(n 2 ) 



(3.10) 



by means of (|3.l| ), with to is given implicitly in terms of z\ and t\ by (3.1), and 
with Kv^(t) defined by an equation corresponding to ( |3.6| ) but for 7 in place of 7: 
clearly one has nvS.it) = KV 1 (t) + 0(k 2 ). For nearby sources (3.10) gives 



St Q 



KVjiti) 

H 2 (h) 



zi + 0{zl) + 0(k 2 ) for zx < 1 



(3.11) 



by means of (|3.2|). 

One may define a spatial proper distance of the source from the observer, with 
respect to the metric g, by 



so ■= c 1^ {S^ — ^-Y^dX = dh t ) 
Jx=x„ d\ 



(3.12) 



where (3.3) has been employed on the right. For nearby sources (3.2) and (3.11) 
combine to give 



ti — tn 



Zl . (i _^M) +0( ^2 )+0(k2) forZl<<1 



H^) Hih) 



(3.13) 



so, for such sources, the redshift Z\ of the source is expressible as a function of its 
spatial proper distance so according to 



, l(so ) = Bb>0 (1 + ^5M) + O(sl) + 0(« 2 ) for Zl « 1 

mi) 



(3.14) 



Notice that in the case Vj{t\) = O(k) equation (3.14) agrees, to first order in Z\ 
and k, with the corresponding formula for the unperturbed FRW space-time. 
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4. Luminosity distance and apparent magnitude 

Suppose now that, in the perturbed space-time, the source with world line {x = 
Xo} radiates uniformly in all directions with power P. 

Let k a := ^2— be the tangent to 7 and let m be a complex vector at (to,Xo) 
satisfying mm = 1, mk = mm = mu = 0. In order to maintain these conditions 
along 7 one requires that m be propagated along 7 according to 

V k m=-(^)k. (4.1) 
v uk ' 

The congruence of all the null geodesies emanating from (to, xo) gives rise to a family 
of Jacobi fields along 7 with expansion g := k a j,m a m b and shear <j := k a .j 1 m a m b with 
respect to the 2-frame {m, rh} and the afhne parameter A. Since the congruence is 
irrotational, the imaginary part of g is zero. The standard propagation equations 
for g and ? therefore reduce to 

± e = -?-<z-e (4.2) 

jLs=-2qg-9 (4.3) 
aX 

where 

6 := ^R ab k a k b = X -KT ah k a k h (4.4) 
* := Rabcdk a k c m b m d = C abcd k a k c m b m d (4.5) 
are the Ricci and Weyl scalars respectively. (The fact that m is propagated along 



7 according to (4.1), rathe r tha n bei ng p arallclly propagated as is more usual, has 



no effect on the equations (4.2) and (4.5).) 



Consider a narrow beam of light rays, centred upon 7, from the source point 
(io,xo)- Let A A and / be the cross-sectional area and apparent luminosity of the 
beam as determined by comoving observers situated along 7. A standard argument 
from photon conservation gives that the quantity (1 + z) 2 IAA is constant along 7. 
By means of g = | In AA one thus obtains 

a = Hrl/2 ^ = Tx ln(/ " /2(1 + zyl) ■ (4 ' 6) 

In terms of the quantities 

J{X):=r 1/2 {l + z)- 1 (4.7) 

e(A):= j2(a) (J9 ^ (a) (4 - 8) 



equations (4.2) and (|4.3|) become 



^ J(A) + J- 3 (A)«A)e(A) + J(A)6 = (4.9) 

J.(| £W ) + ^(A,* = 0. (4.10) 



CUMULATIVE GRAVITATIONAL LENSING 



10 



Substituting for A in terms of t in (4.9) and (4.10) by means of ([Tf]) one thus 
obtains 



d 2 T , v 
dT^ 



S(t) 



+ KI/ 7 (t) J -|j(t) + £(t)£(i) J~ 3 (t) + KT^(t)J(t) = (4.11) 



£(*) - ( IS + ««*r(*) ) £(*) + J 2 (t)K^(t) = (4.12) 



for 



'7 



(*):=o«T, 



a 6 ' 



S(t) 



~dT~dt 



i ,.\ ^ dx a dx c — u — a 



s 2 (t) 



(3) r 1 k lj Bi k m J m 1 ^- 



(4.13) 
(4.14) 
(4.15) 



The initial conditions for (4.11) and (4.12) are 

£(*o) = 
J(i ) = 



d_ 

dt 



J(t) = C1 



t=tn 



(4.16) 



To first order in k, the solutions to (4.11) and (4.12), subject to the initial 
conditions (4.16), are 



J(t) 



P'Sttoj X =t S{t ' ){1 + + °( K ) 

~t'=t V 2, 



4^c 2 S(t) /- t= * £ 2 (i ,i') 



^kV 7 (^)^' + 0(k 2 ) 



foi 



and 



P S(t )J t ,= t0 S(t>) 



rt' = t 

E(*o,t) := / 5(t')d*' 



(4.17) 
(4.18) 

(4.19) 

(4.20) 



where KT^{t) and Kif>^(t) are defined by equations analogous to ( 4.1 3| ) and (4.14) 
for 7 in place of 7. Clearly one has KT^(t) = «;r 7 (t) + 0(k 2 ) and Ki/}%(t) = Kipj(t) + 
0( K 2 ). 

From Q4.17|) and ( |4.18| ) one obtains 



S(t )S(t) dt f =t S 2 (? ,t') 



£ 2 (i ,i) dA ^*'=*° 
for the shear as a function of t. 



'j^ K ^°(t')dt' + 0{k 2 ) (4.21) 
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Writing to — to + Sto and setting t = t\ in ( 4.17 ) one obtains 



ffi+ii£ {s( ; , (l ,(^) s( ? , fl)) 

F 6 W S(t ) 



+ / o S(t)K(jq(t)dt} + 0(0 (4.22) 



by means of ( 3.10 ) and 4.19 ), where to is given in terms of t\ and Zi by (3T). From 
this one obtains 



P S(h) 



S(ti) J t =i 



}(t)dt 



„ £(t ,t)(/^(i)--^ 



KT^t^dt 



S(t ,*i) 



(1 S(t ,t) 



;(£)(&} + 0(k 2 ) (4.23) 



fc ? 5(0 

by means of ( 4.20| ), ( |4. 19|) and an integration by parts. For nearby sources ( 4.23 ) 
gives 

I- 1,2 izi) = M4^^-'¥ hl )+0{zl) + 0{J) forz 1 «l (4.24) 
V r H{tx) H(h) 

by means of ( 3.13 ). The equation corresponding to ( 4.23j ) for the unperturbed 

o . 

metric g is 



-1/2/. , _ „./ 47r ( 1 + z i) 2 w? 



P S{h) 



S(t ,*i) ■ 



(4.25) 



by 



From (4.21), the total shear of the image at the observation point (fi,xi) is given 
A=Al ' ,t=tl S(t )S(t) 



sd\ 



A = An 



"* ^ 2 ( to j^l K ^°(t')dt' I dt + 0( K 2 ) . 



/t=t E 2 (? ,i) y^t'=*o 

The luminosity distance of the source is defined by 



d L := 



4nl(zi) 



(4.26) 



(4.27) 



which, for nearby sources, combines with (4.24) to give 



d L = 



CZ\ H KV°{tij 



(1 - Zp^iZ) + 0{zl) + 0{k 2 ) for zi « 1 . (4.28) 

Note that in the case ^(ii) = O(n) this agrees, to first order in z\ and k, with 
the corresponding formula for the unperturbed FRW space-time. From (4.28)) and 
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( 3.14 ) one has 

d L = s + 0(s 2 a ) + 0(k 2 ) (4.29) 

which gives the acceptable result that the spatial proper distance and luminosity 
distance agree to first order in sq and n, independent of the value of v^ti). This 

o 

holds for both the perturbed metric g and the unperturbed metric g. 

The apparent magnitude of the source at x = Xo relative to a reference source 
at x = x re f is defined by 

5 If 

m := m ref + - log 10 -j 1 (4.30) 

where I and I re f are the apparent luminosities of the respective sources. The refer- 
ence source is taken to have power P re f and to be at a spatial proper distance s 0ire / 
that is small on a cosmological scale (10 parsecs is conventional). For the unper- 

o 

turbed metric g, a source of apparent luminosity / has an apparent magnitude 

o 

o o 5 1 ref / o -i \ 

m = mref + -z lo gio ~ir ( 4 - 31 ) 
z I 

o 

relative to a source of apparent luminosity I re f. The sources will again be taken to 
have powers P and P re f respectively, with the reference source at the same spatial 
proper distance Sq re f as for the perturbed metric g. Let m re f = m re f. The objective 

sources for g and g are both taken to be at redshift Z\. From ( 4.3C ) and ( 4.31 ) one 
obtains 

I ref (s) I(Zl) 

in the limit so.re/ — ► 0- B y ( [4.29), th e first term in the braces is just \a{P/P re f) to 
first order in k. By ( 4.23| ) and ( 4.25 ) one thus obtains 

lnl ° 5(*i)E(i ,ti) Jt=t « Jt = u > b(t > 

' f 1 S(t ,f)(/o^(t) - ^fe^ ^(t))dt> + 0(k 2 ) (4.33) 



S(*o,*i) 



S(t) 



for the correction to the apparent magnitude-redshift relation. Note that this cor- 
rection depends on the perturbation Sh only through the functions Kv$(t) and 
n8TSf(t). 



For nearby sources ( [1.33| ) gives 

Sm(z 1 ) = - - 5 in K f (tl) +0( Zl ) + 0{ K 2 ) forz«l (4.34) 
lnl ° H^) 

for the apparent magnitude-redshift correction. This evidently vanishes at z% = 
to first order in k in the case f 7 (ii) = O(k). 

It is evident from the analysis in the present and preceding sections that the 
equation 

K^ih) = + O(k 2 ) (4.35) 
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may be interpreted as a condition on the perturbation which, as regards local optical 

o 

properties, ensures that the background FRW metric g is a best fit to the perturbed 
metric g, in the direction of j(ti), at the observation point (ti,Xi). A requirement 
that ( 4.35| ) holds for all null geodesies 7 through (ii,Xi) is, by (|3.6|), equivalent to 



the condition 

dtShijih,*!) = + O(k 2 ) . (4.36) 

In general this represent a physical constraint on the perturbation at (£i,xi) since 
there is, in general, no freedom to choose a new time slicing such that the 3-metric 
h is intrinsic to the level surfaces of the new time. However, for the specific class of 
perturbations to be introduced in §||, there is just such a freedom, at least to first 
order in k, which may be exploited to ensure that dOe| ) does hold. 



5. Newtonian perturbations 

In order to describe the matter distribution of the cosmos, physical considerations 
suggest that one seeks a perturbed 3-metric h such that the corresponding space- 
time metric g has an energy tensor of the form 

T ab = (p + p)u a u b +pg ab (5.1) 

describing a perfect fluid with velocity u a :— t' ,a , density p and pressure p. It is 
conventional to define a dimensionless density parameter by 

n=-$. (5.2) 

Note that all of p, H and fl depend on all four coordinates. With regard to the 



pressure, since the form of (2.13) implies that u a is geodesic, the conservation 
equation for (5.1) implies h a b p-b — and hence h a {Sp)-b = under the assumption 
that p + p is nowhere zero. This leads one to consider the particular case Sp = 
which describes perturbations arising from the addition or removal of comoving 
dust. 

The full non-linear problem presents formidable difficulties, although one does 
have by (5.1) and ( 2.25| ) that Shij satisfies the simple equation 

h jk (d t 5h k[i ) m = . (5.3) 

In order to make progress, only a linear approximation to a solution will be sought. 
Specifically the problem is to obtain Sh^, regarded as a power series in k, vanishing 
to zeroth order, such that the energy tensor of g a b has the form 

T ab = (p + p)u a u b + pg ab + O(k) . (5.4) 

o 

Consider a perturbation of the 3-metric hij of the form 

Shij = F(t)K$ K + G{t)K§hij + 0(k 2 ) (5.5) 
for functions F(t), G(t) and a i-independent scalar field ^(x). The form of (|5.5|) 



corresponds to the synchronous gauge perturbations considered by Mukhanov et 
al. (1992, p. 216). It will be convenient, although not necessary, to regard F(t) 
and G{t) as power series in n even though only the zeroth order terms will be of 
significance. Note that although Shij enters into the right side of ( |5.5| ) through the 
term &\\ij, this is not significant in the linear approximation since $>\\ij and $1]^ 
agree to zeroth order in k. 
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From first principles one has 

5 (3) i?.y - 2k8h ik = --G0){K$uy + (k (3) A$ + 4kK$)hij} + 0{k 2 ) (5.6) 



which by ( [2.28 ) gives 

5^R= -2G(t)(K (3) A$ + 3fcK$) + 0(k 2 ) 
Substitution of into (|2.30D , ( |2.3l| ) and yields 



(5.7) 



KT tt -K P +{ s2{t) s2{t) )K A<t> + 3{ s{t) s2{t) ) K <P + U{n ) 



kT« = kTu = (kF(t) - G(t))K#|n + 0(k 2 ) 



(5.8) 
(5.9) 



{S 2 (i)^ + ~(G(t) - J-(S a (t)F(t)) )^A<S> 



(kG(t) 
1 , f 



1 



S(t) 



{S 3 (t)F(t)j -G(t))K$ K + 0(K 2 ) 



(5.10) 



by means of Q, (Fill) and Q . 

For KTi t and the trace- free part of «Ty to vanish to first order in n, in accordance 
with (5^), it suffices to require that F(t) and G{t) satisfy 



From these one obtains 



G(t) = ^ f) (S 3 (t)F(t)) +0(k) 
G(t) = kF(t) + 0(n) . 

(S(t)(S 2 (t)F(t))) =0 + O(n) 



by means of (2J), and hence 

S(t)(S 2 (t)F(t)) =C + 0{k) 



(5.11) 
(5.12) 

(5.13) 

(5.14) 



for some constant C. 

If C were chosen to vanish to zeroth order in k then, by (5.14) and ( |5.11 ), both 
(S 2 (t)F(t)) and G(t) — S(t)S(t)F(t) would vanish to zeroth order in n and so, 
as will be evident from (|6.8| ) in g would be an FRW metric to first order in n 
irrespective of the function $(x). To avoid this uninteresting case C will be chosen 
to be non-zero to zeroth order in n. One may then normalise F(t) to give 

(5.15) 

(5.16) 
(5.17) 



C= 1 



Equations ( |5.1l| ), (|5.14|) and (|5.15| ) give 

(S 2 (t)F(t)) 



1 



S(t) 
1 



O(k) 



G(t)-S(t)S(t)F(t) = — + 0( K ) . 
These are equivalent to equations (5.11), ( |5.12 ) and ( |5.15 ) by virtue of fl2.9|). 
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By means of (5.16) and (5.17), equations (5.8), ( |5.9| ) and ( |5.1C ) combine to give 
that T ab has the required form (5.4) for 



From ( |5.18 ) one has 



5p = + 0{k) . 



d t (S 3 (t)Sp) = + O(k) 



as one w ould e xpect since the p erturbation has a dust equation of state. 
From (|J), (HI) and dUl) one has 

Wdtdhij = -S 3 {t)F{t)n8p + 0{n 2 ) 
whereby ( ^21] ), ( jgjg ) and give 

6H = -ls 3 (t)F(t)KSp + 0(k 2 ) 
6 

Sq = l{S(t)(S 2 (t)F(t)) + 2S 3 (t)S(t)F(t)(l + q(t))}6Q + 0{k 2 ) 
S (3) i? = 2S 3 {t)G{t)n8p + 0{k 2 ) 



(5.18) 
(5.19) 

(5.20) 

(5.21) 

(5.22) 

(5.23) 
(5.24) 



with the help of (|5_2j) and (|5_17j). 

It is straightforward to check, by m ea ns of ( |2.42 ) and (2.11), that the perturbed 
space-time metric g defined by ( [2.13 ), ( [2.14 ) and ( |5.5| ) is silent to first order in n 
in the sense of 



+ O(k 2 ) 



(5.25) 



The electric part of the Weyl tensor of g is, by fl2.39[ ), (|5.6|), (|5_7j), ( |5.16| ) and fl5.17| ), 

given by 



-L-^ ]]l3 -h ljK ^A^ + 0( K 2 ) . 



(5.26) 



by 



By (5^4) and the analogue of ( [4.13 ) for 7, the function Sr^t) in ( 4.33| ) is given 



6T°Jt) = h P (Zf(t))+0(K) 



(5.27) 



An equation for the function v%(t) in ( 4.33] ) will be given in §|| By ( |5.25 ), ( |5.26 ) 

ction ip^(t) in equ 

$ll y m i TO i + 0{n) 



and the analogue of ( 4.15 ) for 7 the function ip^(t) in equation ( 4.26 ) is given by 

1 



S(t) 



(5.28) 



where m % is the analogue of m 1 for 7. 

Space-times satisfying the condition B ao = have been termed 'silent' by some 
authors and 'Newtonian-like' by others. However it is known that the vanishing of 
B a b can be conserved in time only in specialised cases (Maartens et al. 1998). In 
any second or higher order study of perturbations of FRW cosmologies one would 
therefore not expect B a b to vanish to any higher than first order in k. 
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6. The gauge condition 

The functions F(t) and G(t), which enter explicitly into the perturbed metric 
g, are not determined uniquely by ( 5. 16| ) and ( 5.17 ) since these equations allow 
two freely specifiable constants of integration. In alternative terminology, there is 
a gauge freedom 

F(i) - F(t) +A + A ^ (6.1) 

G(t)-G(t)-4>!H (6 ' 2) 

where Aq and A are real constants. In order to understand the meaning of this 
freedom, consider the metric g, as given by ( |2T3|) J2l3| ) and (|J), expressed with 
respect to the coordinate system employed in ([2.2[): 



g = -dt 2 + S 2 (t) (l (1 + G(t)«*(x)) ^(dx 1 ) 2 

+ S 2 (t)F(t)K$(x) Uj dx i dx j + 0{n 2 ) . (6.3) 
In terms of new coordinates (t, x) defined by 

t := t + iF*(t)/s$(x) (6.4) 

1 / At 2 \ 2 

x' ':= + -F(t) (l + — J «6Si*(x), (6.5) 

for an as yet unspecified function i* 1 *^), the metric g assumes the form 
g = -(1 - F*(i>$(x))dP + (F*(i) - S 2 (t)F{t))nd i &(i) dx'dt 

+ ^ (1 + (G(t) - ||V(*))^(x)) ][>*T + 0(* 2 ) . (6.6) 

With the choice 

F*(t) - S 2 (t)F(t) (6.7) 
the metric g assumes the diagonal form 

g = -(l-(S 2 (t)F(t)) K$(x))dP 
S 2 (i) 



(1 ' 4 



r(l + (G(t) - 5(t)S(t)F(t))/c*(x))^(di i ) 2 +0(k 2 ) . (6.8) 



No assumptions about the ^-dependence of F(t) and G(t) have yet been em- 
ployed. It is evident that if (S 2 (t)F(t)j and G(t) - S(t)S(t)F(t) were both to 
vanish to zeroth order in k then, as quoted in §^J, g would be an FRW metric to 
first order in k. 

For F(t) and G(t) satisfying ( 5.16 ) and ( 5.17| ), equation ( |6.8| ) reduces to 
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which is the metric of Newman & McVittie (1982). Note that the curves {x 
const.} have unit tangent 



k$(x)\ - 



0(k 2 ) (6.10) 
They are therefore not the world lines 



and so are non-geodesic to first order in k 
of freely falling observers. 

The absence of F(t) and G(t) in ( |6.9D shows that the isometry class of g is 
unaffected by the choice of integration constants for ( 5.16| ) and (5.17). Indeed the 
gauge transformation (|6.1[), (|6.2|) is induced by the coordinate transformation 



t+^ K $(x) 



A 



df 



kr 2 



(6.11) 



(6.12) 



which preserves the form of the metric (16 



S\f)' 

) for F(t), G(t) transforming according to 



( |Q| ) and (6.2). Moreover the transformation ( |6. 11 ) of t preserves the form of ( p.5[ ) 
on the level surfaces of t, fo r a given <f> o n the space-time manif old, fo r F (t), G(t) 
transforming according to (ST) and (5.2). The gauge freedom ( |6.l[) , ( |6.2|) is thus 
to be int erp reted as a freedom to choos e the ti me fu nction t without violating the 
form of (|5.5|) . The coordinate freedom ( 6.11 ), ( 6.12 ) corresponds to that identified 
by Mukhanov et al. (1992, p.216). 



Although the integration constants for ( 5.16 ) and ( |5.17 ) do not affect the isom- 
etry class of g to first order in n, they nonetheless carry physical significance since 
they help determine the physically significant coordinate t. In order to make phys- 
ically appropriate choices of these constants, consider the function v 1 



along a light ray 7 through the observation point (fx, Xi). By means of ([ 
and the null geodesic equations (|3.3|) and (fTj) one obtains 



of (|3.6) 



^(t) = F(t){S(t)j(s(t)j$( 7 (t)) 



•$( 7 (t))}+0(/s) 
By d5.5| ) and ( |5.12j ) the gauge condition ( 4.36| ) becomes 







This holds for a unique value of Aq in (6.1), (6.2). 
Equations ( |5.16j ), ( |5.17j ) and ( |6.14| ) give 

r t,=tl dt' 



G(t) = 



~S 2 (t)J t , =t 
1 S(t) 



t , =t 1 



O(k) 



(6.13) 
(6.14) 

(6.15) 
(6.16) 



Thus, to zeroth order in k, G(t) is now uniquely specified whilst F(t) is determined 
only up to an arbitrary additive constant. By ( |6.11 ) and ( 6.12| ) this remaining 
gauge freedom corresponds to a coordinate transformation of the form 

t^t (6.17) 

(6.18) 



A 
~2 



kf_ 
4 
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where A is an arbitrary constant. Since this preserves the level surfaces of t it 
has no significance for the physical properties of the perturbed cosmology. There is 
therefore no need to specify the remaining integration constant of (5.16) and (5.17). 

By (4.34) the gauge fixing condition (3.14) may be interpreted as a necessary and 
sufficient condition that there is a vanishing correction to the apparent magnitude- 
redshift relation, in all directions, to zeroth order in the redshift and first order 
in k. By (4.28) another interpretation is that there is a vanishing correction to 
the luminosity distancc-rcdshift relation, in all directions, to zeroth order in the 
redshift and first order in k. 

An immediate consequence of the gauge fixing condition (6.14) is, by (5.22), that 
the Hubble expansion H of the perturbed cosmology satisfies 

fT(ti,x) = H(h) + 0(k 2 ) 



(6.19) 



and so is unperturbed to first order in k at all points of the surface {t = t±}. By 
the use of (5.14) and (6.14) in (5.23) one also has 



Sq 



~sn + o{k 2 ) 



(6.20) 



The gauge condition ( 6.14 ) thus ensures that the perturbation (5.5) affects the 
deceleration parameter but not the Hubble parameter at the observation point 

(tl,Xl). 

Subject to the gauge fixing condition ( |6.14| ) the apparent magnitude-redshift 
r elati on is given by ( 4.33 ), with v%(t) give n by ( 6.13 ) and ( 6.15 ) and Sr^(t) by 
( 5.27). Substitution for these quantities in (4.35) yields 



8m{z\) 



5 f S 4 (t )F(t ) d_ 



S 2 (to) 



11110 L 25(ti)£(t ,ti) dl 



«*(7(*)) 



S 2 (t ) 



S(ti)£(t 0) ti 



S(ii)£(t ,ti) 

rt=t 



1 



S(to, t\ 
t=tl S(to,t) 

o 

t=t 



(1 



S(t) 



t=t 



2S(t) 



S 2 (t) 
K,5p(j(t))dt + - 



E(t ,*))/s$(7(*))<ft 



:tl S(fi_ 
t=i S 2 (t) 

«$(7(*i)) 
2S{h) 

S 3 (t )F(t ) 



2E(t ,ti 



«*(7(*o)) 



S(t ,*i) Jt=t ° 



H X 2 (t ,t) 
2S(t) 



KSp{j(t))dt\ + 0(k 2 ) 
(6.21) 



by successive integration by parts and the use of (5.11), (5.12), (5.16), (5.17) and 



(3.14). Again to is given in terms of z\ and ti by means of ( |3.l| ). Equation (6.21) 
is the correction to the apparent magnitude redshift relation for Newtonian pertur- 
bations of FRW cosmologies subject to the gauge fixing condition (6.14). 
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7. The averaging procedure 

Suppose $ is smooth on each level surface of t. Let Q> t be a 3-domain in one 
such surface. By definition, the mean mass density perturbation on is 



f 5p dv 
f dv 

J S>t 



(7.1) 



where dv is the elemental 3- volume on @t with respect to h. By means of ( |7.l|) and 
( Ljj ) one obtains 



0(«) 



(7.2) 



by means of the divergence theorem, where n is the unit outward pointing normal 

o 

to 3t t at dSlt with respect to h in the surface {t — const.}, and where 



^ 

f dv 



(7.3) 



is the mean value of $ on & t - Note that ( < f) @ is independent of t. It will be 
assumed that the perturbed matter distribution is sufficiently uniform on the large 
scale that, in the limit 9-t — > oo, the surface integral term in ( |7.2| ) tends to zero, 
whilst (bp) ^ and (^)^ tend to limits {Sp) and ($) respectively. Under these 
conditions one obtains 



3fc 



;«> + o(«) 



(7.4) 



In terms of the dimensionless density parameter of (5.2) this gives 



K ($p) t 

3H 2 (t) 



+ 0( K 2 ) = (l + q(t)y 



5(t) 



(7.5) 



by means of (2.9). Equation ( J7.4[ ) and the right side of (7.5) must be considered 
invalid for k — since in that case ($ ) must be infinite to give a finite (5p ) . 
In order to obtain the mean correction to the apparent magnitude-redshift rela- 



tion, one replaces the function 3>(~f(t)) on the right of (6.21) with its mean value 
($ ), and o $(7(f)) by its expectation value of zero. In practice, since all cos- 

mological sources lie within gravitational potential wells, the o $(7(4)) term 



in (3.21) will always give a positive contribution to the apparent magnitude, but 
this will be so small as to be negligible for present purposes. By means of (5.20), 
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(7.4) and (p!q), and for i given in terms of t\ and z\ by (pjj), one thus obtains 



(8m(z\) 

/ 

S*(t)" 



S(ti)E(t ,ti) 

O 

' ^ ~ * E(t ,t)F(t)dt i ^j^ dt } + 0{k 1 



£(to,*i) Jt=to S(*o,*i) Jt =*° 



(7.6) 



for the mean correction to the apparent magnitude-redshift relation for Newtonian 
perturbations of FRW cosmologies. This is the fundamental equation of the paper. 



It could alternatively have been derived directly from ( 4.33 ) w ith th e help of the 
observation that the first term in the braces on the right of ( 6.13 ) gives a zero 
contribution to (8m(zi) ). For nearby sources (|7.6|) gives 



<M*0 > - -TTTTT^ > tl *i + 0(4) + 0(n 2 ) (7.7) 



4 In 10 

by means of ( 3.1 3| ) and (7.5). 



In the particular case of a zero mean density perturbation, (8p) — 0, equation 
( [7.6| ) clearly gives that there is a zero mean correction to the apparent magnitude- 
redshift relation to first order in k. In the case of a non-zero mean density pertur- 
bation (Sp) 0, one may compare either with the background metric g or with 
that of the homogenised cosmology determined from g by the constant potential 
($). (See §||a.) In the latter case the mean perturbation of the potential is zero, 
so the mean density perturbation is also zero. Hence the correction to the apparent 
magnitude-redshift relation is zero, to first order in k. 

It is worthwhile to consider the case (dp) > in more detail. Of the four terms 



in the braces on the right of (7.6), the first three are positive. And though the 
fourth term is negative, it is evidently dominated by the second. The sum of the 
terms in the braces is therefore positive. This shows that an object at a given 
redshift appears brighter than in the reference FRW model. Moreover an object at 
a given redshift in a high density universe appears brighter than in a low density 
universe. These effects can be regarded as the gravitational lensing of the universe 
as a whole. 

A fruitful approach to clumped matter perturbations of FRW cosmologies, in- 
troduced by Dyer & Roeder (1973), is to assume that of all the matter present, a 
proportion a is uniformly distributed and pressure-free, while the remaining pro- 
portion 1 — a is gravitationally bound into clumps. In order to compute the lensing 
of light beams that remain far from all clumps it suffices to take into considera- 
tion the effect of only the uniformly distributed matter. In particular, shear can 
be neglected. The angular diameter distance (Schneider et al. 1992, eq. (3.66)) 
determined from such clump-avoiding light beams is often called the Dyer-Roeder 
distance. For light beams of large angular diameter, one expects the cumulative 
lensing effect to approach that of the homogenised FRW cosmology. The transition 
between these two regimes has been studied by Linder (1998) and found to occur 
typically between 1 and 10 arcseconds. 
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The Dyer-Roeder method may be implemented in the present framework as 

o 

follows. The background metric g is taken to be the metric of the homogenised 
FRW model. The perturbation <5h then corresponds to the removal of a proportion 
a of the matter of this model, followed by the addition of the same quantity of 
clumped matter. The apparent magnitude-redshift relation for narrow, clump- 
avoiding light beams is then given by (7.6), with (Sfl ) replaced by — a(fl ) . By 
means of (2.8) this yields 

(Sm(zi) ) = 



2 In 10 



S 3 (t 1 )H 2 (t 1 )a(SQ[ 



S 2 (t ) 



F(t)dt + 3 



3 



s 4 (t) 



dt. 



T.(t ,t)F(t)dt 



E(t ,ti) J *='° 



S(*o,*i) Jt=t0 



s 2 (t ,t) 



dt 



0( K 2 



(7.8) 

Since the expression in the braces is positive for Z\ > 0, this correction is positive 
for Z\ > 0. Sources viewed at a given redshift along clump- avoiding light beams 
thus appear dimmer than for an all-sky average. This is as one would expect since 
clump-avoiding light beams are less focussed. For larger angular scales, the mean 
correction is no longer given by (7.8) and should be expected to tend to zero since 
the mean density perturbation is zero. Thus (7.8) may be interpreted as the mean 



correction to the apparent magnitude-redshift relation for narrow, clump- avoiding 
light beams relative to that for wide angle beams. 

The standard implementation of the Dyer-Roeder ansatz (e.g. Schneider et al. 
1992, p. 138 et seq.) gives the Dyer-Roeder distance as a solution to the Dyer- 
Roeder equation which describes light propagation through a uniformly under- 
dense region of space-time. From the Dyer-Roeder distance one can obtain an 
apparent magnitude-redshift relation. On the other hand, an apparent magnitude- 
redshift relation for light propagation through a uniformly underdense region of 

(see 



space-time is also described by ([77 

s7 



for an appropriately valued constant ($ 
Nonetheless one cannot expect these two relations to agree unless similar 
gauge fixing conditions are applied in each case (see (8.3). 



8. Examples 

Example 8.1 (Uniform density perturbations). In the special case 

$ = const. 

the perturbation fl5.5| ) reduces to 

Shij = KG(t)<$>hij + 0(k 2 ) . 



(8.2) 



To first order in n this is equivalent to leaving the 3-metric h in (2.1) fixed and 
perturbing the scal e fac tor S(t) according to S(t) — > (1 + ^G(t)n$)S(t). The 
coordinate freedom ( 6.11 ), ( 3.12 ) r edu ces to the freedom to change t by an additive 
constant. One sees directly from ( |8.2|) that the gauge condition ( 4.36 ) is satisfied 
iff G(t) satisfies 

G{ti) = . (8.3) 
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From (|J), Q and (pl| ) one has that the perturbation is pressure-free to 
first order in k, and so of a form as discussed in §B if G(t) satisfies 



kG(t) = —(S 3 (t)G(t)) +0(k) 



(8.4) 



If G(t) were to vanish a t t = t\ then, by ( |8.3| ) and ( p.4[ ), G(t) would vanish for all 
t. For consistency with ( 6.1 6| ) one may choose 

G (^l) = 777TT • (8-5) 



By ( |5.18| ) and ( |5.19| ) one has 



5p = - 



3A- 



$ + O(k) 



s 3 (t) 

Sp = + O(k) . 

By (|7.5|) the perturbation of the dimensionless density parameter is given by 



(8.6) 
(8.7) 



Example 8.2 (Point particle perturbations). A perturbation by the introduction 
of a family of comoving point particles with world lines {x = x a }, a = 1, 2, . . . and 
masses m Q > is described by a density perturbation of the form 



(8.9) 



where ^S Xa (x) is the Dirac distribution with respect to h, centred on x Q , on each 
level surface of t. The par ticles shall represent the galaxies. 

In order to satisfy (|8.9| ) and ( 5.18j ) one seeks a potential $(x) satisfying 

-( (3) A$(x) + 3fc$(x)) = 5^ m Q (3) L Q (x) . (8.10) 

a 

Since this equation is linear one may decompose ^(x) as a sum 

$(x)=^$ Q (x) . (8.11) 

a 

In order to consider a typical summand ( f ) Q (x) it is convenient to express the 3- 

o 

metric h in the form 

' ±(duj 2 a + sin 2 uj a dn 2 a ) if fc > 

„2 ^o2 if = 

^ {duj' 2 a + sinh 2 w Q ) if k < 
where <if2 2 is the 2-sphere metric and the radial coordinate u> a is defined by 

if k > 



h = < dr 2 a + r 2 a dVL 2 a 



12) 



w Q (r Q ) 



2 tan" 1 (^r c 



2 tanh 



1 (^r a ) 



if fc < 



.13) 



with the origin uo a = being the point x a of (S.9). The range of tu a is < uj a < tt 
if k > and < w Q < oo if k < 0. 
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The general radial solution to (S.10), as found by Newman & McVittie (1982), 



is 



*a(x) 



3Ff 



y/k 



-(m a cos 2ui a + C a sin 2uj a ) 



iirr a 



if k > 



if k = 



(8.14) 



-i — —■ — r (m a cosh2w Q + C a sinh2cLi Q ) if k < 

k47rsmna; Q v ; 



where C a is an arbitrary constant. 

The case k > would appear to be the simplest insofar as the level surfaces 
of t are compact, so physical plausibility demands that there are at most finitely 
many particles present. However it is evident from (8.14) in this case that for each 



particle there is a complementary particle of equal mass located at the antipodal 
point in each surface {t = const.}. This bizarre doppelganger phenomenon leads 
one to question whether the k > solution is realistic after all. It is unclear whether 
the problem is an artifact of the symmetry of the level surfaces of t or of the linear 
approximation. There may even be a deeper issue here concerning the constraint 
components of the Einstein equations (D'Eath 1976). 

For k > an integration of (|S.10) over a level surface of t with respect to the 



volume element dv associated with h yields 

-3k / $dv 

J t— const. 



.15) 



Setting ' 3 'y :— J. . dv one thus obtains 



S 3 {t) {3} v 



3k 



for 



Note that 



(4>) = 



f, , §dv 

J t— const. 



dv 



(8.16) 



.17) 



.16) agrees with (7.4) even though <f> is not smooth in the present case. 
In the case k < case the surfaces {t = const.} have infinite volume, so infinitely 
many particles (or none) are needed in order to achieve a distribution that is uniform 
on the large scale. For each a one must choose C a = — m a in order that <fr a decays to 
zero at infinity. Each < I > Q (x) is then integrable on the level surfaces {t — const.} and 
a simple reciprocity argument indicates that for a sufficiently uniform distribution 
of particles ^(x) — J2 a should converge everywhere other than on the world 

lines of the particles. Let be a compact 3-domain in a level surface of t. An 
integration of ( 8.10| ) over S> t with respect to dv yields 



3k $di = Y^ '' 



(8.18) 



try, 



where n is the outward unit normal at d&>t , da is the area element on d£>t induced 

o 

by h, and where, in the first term on the right, the sum is carried over all a such 
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that 



Setting V(@t) ■— J@ dv one obtains 



3k 



J as> t " 



(8.19) 



where 



(Sp) 



53 ( t ) &)v(9 t ) 
is the mean density perturbation on £} t and 



(8.20) 



(3) 



.21) 



is the mean value of $ on 9% ■ Reasonable uniformity conditions on the distribution 
of the particles should ensure that (Sp ) and ($ ) tend to limits (Sp ) and ($ ) 
respectively for arbitrarily lar ge ffi . Such conditions should also ensure that the 
second term on the right of ( 8.19 ) tends to zero for large &t- One then obtains 
(7.4) for k < 0, again even though $ is not smo oth. 

In the k — case $ Q (r Q ), as given by (|S.14 ), decays to zero at infinity only if 
C a = 0, and then only as l/r a . For a uniform distribution of particles the potential 
$ = ~^2 a $ Q would be infinite everywhere, so the theory breaks down in this case. 

For k 7^ the perturbation of the space-time metric corresponding to the in- 
troduction of the uniform distribution of comoving point particles is described by 

(foj), q, _ 



11) and (8.14), and the mean correction to the apparent magnitude- 



redshift relation, relative to the background metric, is given by (7.6) 



Example 8.3 (Swiss cheese model). This model, proposed by Einstein & Strauss 
(1945,1946), is an exact C 1_ solution to the Einstein equations consisting of a 
pressure-free FRW model in which spherical regions are replaced by spherical pieces 
of Schwarzschild geometry. The idea is that each of these spherical 'holes' represents 
the condensation of dust into a star, represented by the singularity at the centre. 
For present purposes, the singularities will be considered to represent the galaxies 
and the intervening dust, the 'cheese', will represent the intergalactic medium. The 
model will be considered here in terms of the perturbative formalism of §||. 

Let p(t) be the density of the intergalactic dust. In order to describe a spherically 
symmetric hole with a central point particle of mass m and coordinate radius u> = Co 
if k =/= 0, or r = f if k — 0, one may seek a radial potential function $(x) which 
satisfies 

-( (3) A$(x) + 3fc$(x)) = -S 3 (t 1 )p(t 1 ) + m (3) 5 (x) (8.22) 



in the hole and matches in a C 2 manner to a constant potential ^(x) = <t outside 
the hole. The first term on the right of (8.22) ensures that the hole is a vacuum. 
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.23) 



The general radial solution to ( 8.22| ) has the form 

' Vkm cos 2uj , C sin 2uj , ff 3 (ti)p(fi) 
47T sinu B7F sin a; 3fc 

$(x)=<!^ F + ^+l5 3 (t 1 )p(t 1 )r 2 

f^km cosh 2a; , C sinh 2u , g^tQpfti) 
47r sinh 2w B7F sinh w 3fc 

within the hole, where C is a constant. 

In order for the solution 

tion <£>(x) = $ outside the hole, the radial derivative of $(x) must van ish at the 

boundary. By means of the divergence theorem, an integration of ( 8.22 ) over the 

hole therefore yields 



if k > 
if fc = 

if k < 



.23) to join in a C manner to the constant solu- 



m = S 3 (t 1 )p(t 1 )^V hole -3k 



(8.24) 



hole 



for all values of k, where ^Vhole /1 . d» is the volume of the hole with respect 

o 

to h. Thus the mean matter density of the h ole is prec isely p{t\) for k = and 
p(h) + 0{Cj 2 ) for k ^ 0. Substitution of (||J|) into dHJ) yields 



C 



-2Vfcm ( ;° S "^f COS o 3 ^ 
(sinw — g sinow) 



-fcm 



( 5 cosh 3u) — cosh u)) 
( 5 sinh 3d) — sinh a)) 



if fc > 
if fc < 



.25) 



whilst in the case k = flS.24 ) gives 



yp(ti)S 3 (ti)r 3 



if fc = 



For k ^ the continuity of $ at the boundary of the hole gives 



S 3 (h)p(ti) ^fkm 



by means of 



~3F 

3fc 



1 



(sin Q) — I sin 3d>) 



—km 
"3^ 



5 sinh 3a) — sinh u>) 



if fc > 
if fc < 



723| ) and (8.25), whilst in the case fc = one obtains 
3m C 



(8.26) 



.27) 



(8.28) 



$ = r + — if fc = 

87rf 47r 

by means of Q8.23p and ( |8.24| ). 

Suppose now that there are many holes, each labelled by an index a. Let 3>t 
be a compact 3-domain in a level surface of t such that dS^t intersects none of the 
holes. The mean value of $ on @ t is 



9>t 



(•3) 



(3) 



V(^ t ) 



(8.29) 



o o 

where ^V(^t) ■— f@ dv is the volume of &t with respect to h and the sum is 
carried out over all a such that the a th hole is contained in @f By means of (3.24) 
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and (3.27) in the cases fc ^ 0, and by means of 
k = 0, one has 



where 



M, 



2>t 



— 



1 



1 



(3) 



1 



l v TO « / 

3 Z^a ~~fcl 



_J3j), ( J8.26D and 
2(d) Q — o sin2w Q ) 



28 ) in the case 



(8.30) 



(sinwc 

2 



■ sin3cD Q ) 



1} 



(•'!) 



3 fc i 1 



2(i sinh2w a — w a ) 
(i sinh3a) Q — sinhd) Q ) 



if fc > 
if fc = 

} if fc < 



.31) 



It will be assumed that the distribution of holes is sufficiently uniform that M@ t 
and ($ )g tend to limits M and ($ ) respectively as Q>t becomes arbitrarily large. 

In order to obtain (<& ) = one must have 



By (7.4) one then has 



$ + M = 



(*p) t = + O(«) 



.32) 



(8.33) 



In the case fc = the quantities rh a and f a are, for each a, related by an equation 
of the form ( [3.24 ) , whereby one has 

2/3 z a mr 



Mm = ( 



3tt 



if fc = 



.34) 



For fc 7^ the quantities rh a and f a are, for each a, related by an equation of 



the form 



27J) which involves the constant <I>. However the holes may be assumed 

and 



sufficiently small that (3.34) is a valid approximation for fc ^= 0. By 
( |3.34 ), the dimensionless density parameter of the intergalactic matter is then to 
be perturbed by an amount 



SQcheese = -(1 + 



kM 



for 



/ 3tt 



M 



2/3 



lim 



if fc f 



5/3 



(8.35) 



(8.36) 



the sum being carried over all a for which the a hole is contained in @ tl , with 
V tl (9 tl ) = S 3 (h) (3) y(^ tl ) the volume of 2> tl with respect to the 3-metric S" 3 (ti)h 
induced by g on {t — t{\. If one regards the fc = case as a limit of the fc 7^ 



cases then 



35) may considered to apply for all fc. 



It is the perturbation (8.35) of the density of the intergalactic matter that distin- 



guishes the present approach to the Swiss cheese model from those of other authors, 
including Dyer & Roeder (1974). Only with such a perturbation will the mean den- 
sity parameter (Q ) be unperturbed. And only then, as discussed in will there 
be a zero mean correction to the apparent magnitude-redshift relation to first order 
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9. Concluding remarks 

The validity of the linear approximation employed in this paper depends upon 
the smallness of the mean dimensionless density perturbation (SQ } associated with 
the matter in the galaxies. Astronomical estimates of this quantity are at present 
inconclusive, although there is a general consensus that one does have (SSI ) <C 1. 
If this is correct then the theory presented should provide a valid description of 
the total gravitational field and cumulative weak lensing effects of the galaxies for 
any given model of the distribution of galactic matter. Depending upon the actual 
value of (Sfl ) and the range of redshifts of interest, there may be a need to carry 
the theory to second or higher order in n. This would at least take into account 
the effect of shear on the apparent brightness of cosmological sources. Caustics 
however cannot be adequately described by any finite order power series analysis. 
For this it would be necessary to consider the full non-linear theory. 

I express my appreciation to Prof. M.A.H. MacCallum for his interest and encouragement, 
and to the Department of Mathematics and Applied Mathematics at the University of 
Cape Town for their hospitality. Funding was provided by the Swedish Natural Science 
Research Council (NFR). 
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